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Abstract. The purpose of this article is to show that a C 1 differential system on 
M. n which admits a set of n — 1 independent C 2 conservation laws defined on an 
open subset $1 C R n , is essentially C 1 equivalent on an open and dense subset 
of £7, with the linear differential system u'i = ui, u' 2 = u%, . . . , u' n = u n . The 
main results are illustrated in the case of two concrete dynamical systems, 
namely the three dimensional Lotka-Volterra system, and respectively the 
Euler equations from the free rigid body dynamics. 
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1. Introduction 

Recently, in [G] it is proved that an integrable C 1 planar differential system is 
roughly speaking (^equivalent to the linear differential system u[ = Ui, u' 2 = u-2- 

The purpose of this article is to generalize this result in the n dimensional 
case for a C 1 differential system that admits a set of n— 1 independent conservation 
laws. In the second section we show that such a system can always be realized as a 
Hamilton-Poisson dynamical system on a full measure open subset of R™ with re- 
spect to a rank 2 Poisson structure. In the third section a new time transformation 
will be explicitly constructed in order to bring the system to a linear differential 
system of the type u' 1 = u\, u' 2 = u<z, ■ ■ ■ , u' n = u n . In the last section we illustrate 
the main results in the case of two concrete dynamical systems, namely the three 
dimensional Lotka-Volterra system, and respectively the Euler equations from the 
free rigid body dynamics. 

For details on Poisson geometry and Hamiltonian dynamics, see, e.g. [1], [2], 
[11], [8], [9], [10], [12]. 
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2. Hamiltonian divergence free vector fields naturally associated to 
integrable systems 

In this section we give a method to construct a Hamilton-Poisson divergence free 
vector field, naturally associated with a given Hamilton-Poisson realization of a n 
dimensional differential system admitting n — 1 independent conservation laws. 

First step in this approach is to construct a Hamilton-Poisson realization of 
a given n dimensional differential system admitting n — 1 independent integrals of 
motion. 

Let us consider a C 1 differential system on R": 

±! = Xi(x!, ...,X n ) 

±2 = X 2 {x 1 , . . . , X n ) ^ 

x n — X n (xi , . . . , x n ) , 

where Xi,X%, . . . ,X n G C 1 (R",1R) are arbitrary real functions. Suppose that 
Ci, . . . , C n -2, Cn-i : Q R" R are n — 1 independent C 2 integrals of mo- 
tion of (2.1) defined on a nonempty open subset fi C R™. 

Since C±, . . . , C„_2, C«-i : ^ C R n — >• R are integrals of motion of the vector 
field X = X\d Xl + ■ ■ ■ + X n d Xn S X(R"), we obtain that for each i <G {1, . . . , n— 1} 

n 

(Va(x),X(a;)) = • = 0, 

i=i 

for every x = . . . ,x„) € £1, where (•, •) stand for the canonical inner product 
on R™, and respectively V stand for the gradient with respect to (•,•). 

Hence, by a standard multilinear algebra argument, the C 1 vector field X is 
given as the C 1 vector field *(VCi A • • • A VC„_i) multiplied by a C 1 real function 
(rescaling function), where * stand for the Hodge star operator for multivector 
fields (see for details e.g. [5]). It may happen that the domain of definition for the 
rescaling function to be a proper subset of £1. In the following we will consider the 
generic case when the rescaling function is defined on an open and dense subset of 
£1. In order to simplify the notations, we will also denote this set by f2. 

Consequently, the vector field X can be realized on the open set 51 C 1™ 
as the Hamilton-Poisson vector field Xh G X(f2) with respect to the Hamiltonian 
function H := C n -i and respectively the Poisson bracket of class C 1 defined by: 

{/;3}i';Ci,...,C„_2 ^i A • • • A dx n = vdCi A . . . dC„- 2 A df A dg, 

where v G C 1 (fi,R) is a given real function (rescaling). For v = 1, the associated 
Poisson bracket in the smooth category, it is exactly the Flaschka-Ratiu bracket. 
For similar Hamilton-Poisson formulations of completely integrable systems see 
also [3], [7]. 

In coordinates, the bracket {f,g}v;C u —,C n - 2 is given by: 

r, , a(c lt ...,c n -2,f,g) 
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Note that {C\, . . . , C n -i\ is a complete set of Casimirs for the Poisson bracket 

{•) ■}v;C 1 ,...,C n -2- 

Recall that the Hamiltonian vector held Xh £ X(STj is acting on an arbitrary 
real function / e C fc (^,R), k > 2 as: 

X H {f) = C , G C 1 ^). 

Hence, the differential system (2.1) can be locally written in O as a Hamilton- 
Poisson dynamical system of the type: 

±i = {xi, fl'}i/ ; a 1 ,...,C„_ 2 
±2 = {x 2 , H} v - t c u ...,c n -* 



or equivalently 



Xl = v 



x 2 = v ■ 



{Xn, H}v t C u ...,C n - 2 , 
d(d,. . . ,C n -2,Xl,H) 

d(xi,...,x n ) 
d(Ci,...,C n - 2 ,x 2 ,H) 

d(xi,. ..,x n ) 
d(C\, . . . , C n -2,x n ,H) 



(2.2) 



d(xi, . ..,x n ) 

Consequently, the components of the vector held X = X\d Xl + • • • + X n d Xn 
which generates the differential system (2.1), are given in il as follows: 

x = v _ d{C 1 ,...,C n -2,x i ,H) 
d(x 1 ,...,x n ) 

for i G {1, . . . , n}. 

Next result gives a method to construct a divergence free vector held out 
of the vector held X. The divergence operator we will use in this approach is 
the divergence associated with the standard Lebesgue measure on R™, namely 
Cxdxi A • • • A dx n = divXcfai A • • • A dx n , where Cx stand for the Lie derivative 
along the vector held X. 

Theorem 2.1. The vector field X := — • X is a divergence free vector field on 



£1 \ Z(v), where Z(v) = {{x\ 



,x n ) e Q I u(xi, ...,x n ) = 0}. 



X„.d~ are 



Proof. Note that the components of the vector field X = X\d Xl + 
given by: 

j£ = d{Ci, . . . 7 C„-2,x t ,H) 
d(xi,. ..,x n ) 

for i £ {1, . . . ,n}. By definition, the vector field X is a Hamilton-Poisson vector 
field with respect to the Flaschka-Ra^iu bracket, and having the same Hamiltonian 
H as the vector field X. 
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Hence, the divergence of X is given by: 

<9(Ci, . . .,C n - 2 ,Xi,H) 



d(xi, ...,x n ) 



. U\X\ , . . . , Xi , . . . , X n ) 



where the notation "aij" means that "x^ is omitted. 

Let us now analyze the general term in the above sum. By using the derivative 
of a determinant wc obtain the following: 



a d(C 1 ,...,C n - 2 ,H) ^ d(d Xi d,...,C n - 2 ,H) | | d(C 1 ,...,d Xi C n - 2 ,H) 
<9(x i ; • ■ • ) Xi , . . . , x n ) d{x \ j ... , ... 5 Xji) 9{xii ■ • ■ i Xi , . . . , x n ) 

| d(C 1 ,...,C n - 2 ,d Xi H) 

(9(xi, . . . , Xj, . . . j x n ) 



Hence, 



div(X) = + . . . + d(Cu.-.,d Xi C n - 2 ,H) 

O^Xi , . . . , X^ , . . . , X n ) 0\X\^ . . . , Xj, . . . , X77, J 

, fl(C lj ...,C n _ a ,Q 1 , 4 F), 



<9(xi , . . . , Xi , . . . , Xn ) 



) 



d(x\ , . . . , x^ , . . . , x n ) 9 (xi 5 ... , Xj , . . . , x n ) 



_i <9(Ci, . . . , C n -2,d Xi H) 

3\X\ , . . . , x^ 5 . . • 5 x n ) 



Next we prove that each of the above sums vanishes. In order to do that, it is 
enough to show that the general sum Sk vanishes, where 



y-^_^i+ n _i d(Ci, . . . , Ck-i,d Xi Ck,Ck+i, • ■ • , C n -2,H) 



d(x\ , . . . , Xi , . . . , Xn) 
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Indeed, we obtain that: 

5 fc = iy +k d 2 Ck ■ ■ • ■ ; Cfci ■ ■ ■ ■> Cn-2, H) 

+E(-ir + - 1 E (-d^- 1 ^^ • ■ • • ft' • • •i cr "- a> f ) ) 

_ y^^_iy+n-i y^(_i^J+ fc g 2 Cfc ' ' ■ ' ^ fc ' ' ' ■ ' ^ n ~ 2 ' 

i=1 d{%l i ■ ■ ■ j Xj , . . . , Xi, ■ . ■ , x n ) 

_l_ y^^_^y+n-i ^_^^+fe-i^2 ^11? : : : ! 1 : : : ' ^ ra ~ 2 ' ^) 

. =1 i=j+l 5(xi, . . . , Xj, . . . , Xi, . . . , X n ) 

= i _ ir+k -i j2 i(-i) i+j d* jXi c k + (-i)^- 1 ^.^] • gf 71 ' • • • ■ • • -f n - 2, f\ 

l<j<i<n °( x li ■ ■ ■ j Xj, ■ ■ ■ , X,, . . . , X„J 



l<j<i<n 

= 0. 



□ 



Remark 2.2. In the symplectic case, each Hamiltonian vector field is divergence 
free with respect to the divergence operator defined by the Liouville volume form. 



3. A unified linear formulation of integrable systems 

In this section we give a unified linear formulation for the n dimensional differential 
systems admitting a set of n — 1 independent conservation laws. The construction 
of this linear formulation makes use explicitly of the Hamilton-Poisson realization 
(2.2) of a differential system of type (2.1). 

Let us now state the main result of this article. The notations and respectively 
the hypothesis are supposed to be the same as in the previous section. 

Theorem 3.1. Let (2.1) be aC 1 differential system having a set ofn—1 independent 
C 2 conservation laws defined on an open subset fl C R™. Assume that there exists 
a C 1 rescaling function v , nonzero on an open and dense subset Qq ofQ, such that 
the system (2.1) admits a Hamilton-Poisson realization of the type (2.2) and the 
Lebesgue measure of the set 

:={x= ( Xl , ...,£„) eO | div(*)(*) • 9 ( 1 /^---»Cn-2 i g) (a;) = j 
I d{Xl,...,X n ) J 

in £1q is zero. 
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Then, the change of variables ((xi, . . . , x n ), t) H> ((wi, . . . , u n ), s) given by 

ui = l/v(xi, ...,x n ) 
u 2 = Ci(xi, . . .,x n )/v{x±, ...,x n ) 

U n -1 = C n -2(X1,. . .,X n )/v(xi, ...,X n ) 

u n = H(x 1 ,. . .,x n )/v(x x ,. .-,x n ) 
ds = - dxv(X)dt, 

in the open and dense subset Qqq '■= £Iq\0 o/Qq, transforms the system (2.2) 
restricted to Qqo into the linear differential system = U\, u' 2 = U2, ■ ■ ■ , u' n = u n , 
where "prime" stand for the derivative with respect to the new time "s". 

Proof. Let us start by recalling that the change of variables (xi_, . . . , x n ) 1 — > 
(u±, . . . , u n ) is of class C 1 on Hqq. More exactly, the change of variables is defined 
on a wider set, namely Qq\£, where 

£ : = i x = (X!, . . . , x n ) e O I — ; r(x) = 

C*\Xl , • ■ • , x n J 

= < x = (xi, . ..,x n ) e il I (l/v) (x) — (x) = 

^ U\X\ , . . . , X n J 

I ( \ ez O I d{\/ V ,C 1 ,...,C n -2,H) f 

= |* = (« 1> ...,x ri )eno| d(^—^ {x)=0 

In order to determine the transformed differential system (2.1) (generated by X) 
through this change of variables, we need the following result. Recall by Theorem 
(2.1) that the vector filed l/v ■ X is a divergence free vector field on fioo- Since 

div(l/V • X) = (V(l/v),X) + l/v- divpf), 

and div(l/v • X) = 0, we obtain that: 

{V{l/v),X) = -l/v-A\v{X). (3.1) 

After the change of variables in fiooi the system (2.1) becomes: 



{l/v • d) = d ■ {V{l/v),X) + l/v ■ (VC U X) = - div(X) • u 2 

^-{1/v ■ C„_ 2 ) = C n _ 2 • (W(l/v), X) + l/v ■ (VC„_ 2 , X) = - div(X) • u n ^ 
dt 



(3.2) 

where we used the relation (3.1) and the fact that Ci, . . . , C n _2, H are conservation 
laws for the vector field X, and consequently 

(VCi , X) = ■ ■ ■ = (VC n „ 2 ,l) = (VH, X) = 0. 



du\ 


d 


~~dt ~ 

du 2 


1 


~dt ~ 


dt 


du n -i 




, dt 




du n 


d 


~dT ~ 


dt 
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Now using the new time transformation ds = — drv(X)dt on i7och the system (3.2) 
becomes: 

<jl\ = Ui 
•L = u 2 



where u' t = — '—, for i G {1,. . .,n\. □ 
1 ds 

Remark 3.2. If the rescaling function v =: v cst . is a constant function, then the 
Lebesgue measure of the set O = = £1 is nonzero in Qq, and hence the as- 
sumptions of the Theorem (3.1) do not hold. In this case, we search for a new C 1 
rescaling function fi defined on an open and dense subset £Iq of 0, such that the 
vector field fi ■ X satisfies the assumptions of the Theorem (3.1). The function /z 
satisfies: 

div(fi ■ X) = (Vju, X) + n- div(X) = (Wfi, X), 

since in the case of a constant function v = v cst . we obtain from Theorem (2.1) 
that 

= div(l/i/crt. • X) = (1/i/ct.) • div(X), 
and hence div(X) = 0. 

Consequently, we have to search for a rescaling function fx such that 

divOu-X) = (Vm,X) 

it is not identically zero in fl . 

The second condition that fi has to satisfy is that the function 

a(l//J,Ci,...,C n -2,ff) 
d(xi, . ..,X n ) 

it is not identically zero in VLq. 

The transformation between the differential system (2.1) and the differential 
system generated by the vector field /j, ■ X is done by using the new time transfor- 
mation dt = [i(x)dt' . More exactly, the differential system (2.1) generated by X, 
namely: 

dx , > 
Tt =X{x) > 

is transformed via the new time transformation dt — fx(x)dt' into the system: 

d ^=n{x)-X{x). (3.3) 

The differential system (3.3) can also be realized as a Hamilton- Poisson dynamical 
system with respect to the Poisson bracket {•, •} f4 .i/ <Mt ;Ci,....C„_ 2 an d respectively the 
same Hamiltonian H as for the Hamilton- Poisson realization (2.2) of the vector 
field X. 
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4. Examples 

In this section we will apply the main result of this article in the case of two 
concrete dynamical systems, namely a 3D Lotka-Volterra system and respectively 
Euler's equations of the free rigid body dynamics. 

Let us start with the Lotka-Volterra system. The 3D Lotka-Volterra system 
we consider (see e.g. [4]), is described by the following differential system: 



dxi 

~~dt 

dx 2 



xx(x 2 + x 3 ) 
x 2 (-xi + x 3 ) 



(4.1) 



dt 

dx 3 

— = x 3 {-xi - x 2 ). 
dt 

If one denote: 

X = [xi(x 2 + x 3 )]d Xl + [x 2 {~xi + x 3 )]d X2 + [x 3 (-xi - x 2 )}d X3 , 

then div(V) = 2(x 3 — x\). 

The system (4.1) admits a Hamilton-Poisson realization of the type (2.2), 
where: 



l/(xi,X 2 ,X 3 ) - 

C{x\,x 2 ,x 3 ) 



X\ + x 2 + x 3 
x 2 {x\ + x 2 + x 3 ) 



XlX 3 

H(xi,x 2 ,x 3 ) = xi + x 2 + x 3 . 

The sets introduced in Theorem (3.1) in the case of the Lotka-Volterra system 
(4.1) are given by: 

Sl = {{x 1 ,x 2 ,x 3 ) £t 3 \xix 3 ^0}, 

= {{xi>X2,x 3 ) e R 3 I x x x 3 ^ 0; x-y + x 2 + x 3 ^ 0}, 
O = {{xi,x 2 ,x 3 ) e R 3 I xi = x 3 ; x x x 3 ^ 0; x\ + x 2 + x 3 ^ 0}, 
^00 = {(xi,x 2 ,x 3 ) e R 3 I xi ^ x 3 ] x\x 3 ^ 0; x\ + x 2 + x 3 ^ 0}. 

Then by Theorem (3.1), the change of variables {{x\,x 2 ,x 3 ),t) 1— >■ {{u\, u 2 , u 3 ), s) 
defined by: 

x\ + x 2 + x 3 



Ul 

"2 

"3 



x 2 (xi + x 2 + x 3 y 

x\x 3 

(xi + x 2 + x 3 y 

x\x\ 



, ds = 2(xi — x 3 )dt, 
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for (xi, X2, x 3 ) € f2oo transforms the Lotka-Volterra system (4.1) into the linear 
differential system: 



s 

U 2 



s 

u 3 

ds 



Ul 



U2 



U 3 . 



Let us now analyze the Euler equations from the free rigid body dynamics 
(see e.g. [8], [9], [10], [11], [12]). Recall that the Euler equations from the free rigid 
body dynamics are given by the differential system: 

dxi _ 1 2 - I3 

dt ~ hh 
dx 2 _ h ~ h 

dt ~ hh 
dx 3 _ 1 1 - I 2 

dt ~ hl 2 

where the nonzero real numbers h, h, h are the components of the inertia tensor. 
In the following we consider the case when h ^ h- 
If one denote: 



-X2X3 



-X1X3 



-X1X2, 



(4.2) 



X 



h 



-x 1 x 2 )d x 



( T T x 2 x 3 )d xi + ( xix 3 )d X2 + ( 

I2I3 -il^3 ^l-'2 

= 0, and hence the assumptions of the Theorem (3.1) do not hold. 



then div(X) 

The system (4.2) admits a Hamilton-Poisson realization of type (2.2), where: 



V(XI,X 2 ,X 3 ) =: V cst .(xi,X2,X3) = 

C(xi,X2,x 3 ) = -(xj +x\ + x\), 



-1, 



H{xi,x 2l x 3 ) 



2 [ h 



— + — ). 

In order to correct the vector field X such that one can apply the Theo- 
rem (3.1), we use the Remark (3.2) and consequently choose a rescaling function 
n(x\, X2, x 3 ) = xi, and the associated new time transformation dt = [i(xi, x 2 , x 3 )dt' 
defined on the open and dense subset of R 3 given by {(xi, X2, x 3 ) & R 3 \ xi ^ 0}. 

This new time transformation, transforms the system (4.2) into the differen- 
tial system: 

dx\ h — h 

~d¥ 

dx 2 

~dV ~ hh 

dx 3 h - h o 

-x\x 2 , 



X\X 2 X 3 



hh 

h h 2 

-x\x 3 



(4.3) 



dt' 



hh 



10 



Razvan M. Tudoran 



Recall that the vector field which generates the differential system (4.3) is jj, ■ X: 



v , h — h \ c\ , fh — h 2 \ o , / h — ^2 ; 

(J.-X = {—j—f — xiX2X 3 )d Xl + {—j—f- — x 1 x 3 )o X2 + {— — — x 1 X2)o X3 . 
HLz H-I3 HI2 



'3 " J l 2 \ o . / J l - J 2 2 
-X 1 X 3 )d 2:2 + (— jt-jt 

One note that the divergence of the vector field /x • X is given by: 

div(^i • X) = 12 / 3 X2X3. 

For I2 ^ I3 we have that div(/i • X) it is not identically zero. 

The system (4.3) admits a Hamilton- Poisson realization of the type (2.2), 
where: 

Z/(xi,X 2 ,X 3 ) = V cs t.(Xl,X2,X 3 ) ■ pt{x\, X 2 ,X 3 ) = -Xi, 

C(x 1 ,x 2 ,x 3 ) = hxl +x 2 2 +xj), 

The sets introduced in Theorem (3.1) in the case of the system (4.3) are given by: 

O = M 3 , 

fi = {(n,i2,i 3 )ei 3 I xi ^o}, 



o = {(xi,x2,x 3 ) e 



z\ / 0; x 2 x 3 = 0}, 



fin 



{(xi,x 2 , x 3 ) e M 3 I XiX 2 X 3 ^ 0}. 



Then by Theorem (3.1), the change of variables ((xi, X2, x 3 ), t') h-> ((u\,U2,u 3 ), s) 
defined by: 

1 



ui 

U 3 

ds 



2xi 



x\ 



2/i 2xi/ 2 2xi/ 3 

3 r T 2 X2X 3 dt', 
I2l 3 



for (xi, X2, x 3 ) G finn transforms the system (4.3) into the linear differential system: 

ds 
du 



U2 

is 

■u 3 

ds 



ds 
di 



U2 
U 3 . 
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